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Thermodynamic properties of quantum sine-Gordon spin chain system KCuGaFg 
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We investigated the thermodynamic properties of the spin-i one-dimensional Heisenberg anti- 
ferromagnet KCuGaFg by measuring the specific heat in magnetic fields. When this compound is 
subjected to a uniform magnetic field H, a transverse staggered magnetic field h is induced in this 
compound owing to the staggered component ottheg tensor and the Dzyaloshinskii-Moriya inter- 
action with an alternating D vector. Consequently, the quantum sine-Gordon (SG) model is an 
effective model of this compound in a uniform magnetic field. In three different field directions, we 
observed a magnetic-field-induced gap, which increases with H. We analyzed experimental results 
using specific heat theory based on quantum SG theory. The thermodynamic property for || c is 
very well described in terms of the elementary excitations characteristic of the quantum SG model, 
while for the other field directions, significant contributions from other excitation modes beyond the 
framework of the quantum SG model were observed. For H \\b, a quantum phase transition between 
gap less and gapped ground states was observed. 
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PACS numbers: 75.10.Jm, 75.10.Pq, 76.30.-v, 76.50.- 

I. INTRODUCTION 

One-dimensional (ID) quantum systems have at- 
tracted considerable attention for a long time because of 
various phenomena resulting from the strong quantum 
fluctuations characteristic of the ID systems. In par- 
ticular, the S=l/2 antiferromagnetic Heisenberg chain 
(AFHC) has been closely studied theoretically. Using the 
Bethe ansatz, the dispersion relation of the lowest exci- 
tation for the S = 1/2 uniform AFHC, which is known as 
the des Cloizeaux and Pearson (dCP) mode, has been cal- 
culated exactly-i-!^. The two-spinon contribution to mag- 
netic excitations has also been calculated.—*^ The energy 
of the dCP mode is different from the result obtained by 
linear spin wave theory^i^ by a factor of tt/2, which was 
verified by inelastic neutron scattering measurements 
Accurate calculations^iiiS and inelastic neutron scatter- 
ing measurementsii showed that low-energy excitations 
in a magnetic field cannot be described even qualitatively 
by linear spin wave theory. The calculated dispersion 
relation demonstrated that gapless excitations occur at 
incommensurate wave vectors q = ± 27rm(i7) = ± go and 
TTibgo in addition to at q = and tt, where m[H) is the 
magnetization per site in the unit of g/^B- 

By means of neutron inelastic scattering and spe- 
cific heat measurements, Dender et alr^ observed an 
unexpected magnetic-field-induced incommensurate gap 
in Cu(C6H5COO)2-3H20, abbreviated to copper ben- 
zoate, which is known as an 5*= 1/2 ID antiferromagnet 
with good one-dimensionality. ^^^^^ This problem was dis- 
cussed by Oshikawa and Afflecki^">i^ on the basis of the 
effective Hamiltonian expressed as 

n = J2 [JS^ ■ - Ql^BHSt - {-lr9^iBhS^] , (1) 

i 

where h is the transverse staggered field induced by 



the external field H. The staggered field on the ith 
site hi originates from the staggered component of the 
g tensor (— l)*gs and the Dzyaloshinskii-Moriya (DM) 
interaction^^ with the alternating Di vector (— 1)'D, and 
is expressed as 



hi 
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where is the uniform g factor Because the gs tensor 
and D vector are anisotropic, the proportionality coeffi- 
cient Cs = h/H depends on the field direction. 

Using the field-theoretical approach, Oshikawa and Af- 
fleck demonstrated that the model given by eq. Ilj can 
be mapped onto the quantum sine-Gordon (SG) model 
with Lagrangian density 



dt 



/iCcos(27ri?0), (3) 



where (j) is the canonical Bose field, (j) is the dual field, 
R is the compactification radius, v is the dimensionless 
spin velocity and C is a coupling constant. The dual 
field 4> corresponds to the angle between the transverse 
component of the spin and the reference direction in a 
plane perpendicular to the external magnetic field. Ow- 
ing to the nonlinear term in eq. all the excitations 
are gapped as illustrated in Fig.[Tl 

The elementary excitations characteristic of the quan- 
tum SG model are the soliton, antisoliton and breather 
excitations. The soliton excitation corresponds to the 
excitation at q = ±qQ and tt ± go i and is classically de- 
scribed as the local rotation of the spin in a plane per- 
pendicular to the magnetic field, and the antisoliton ex- 
citation corresponds to the inverse rotation of the spin. 
The breather corresponds to the bound state of the soli- 
ton and antisoliton as reducing their energies with </) ~ 0. 
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Oshikawa and Affleck-lfi"-l§ showed that the magnetic- 
field-induced gap corresponding to the sohton energy is 
proportional to H'^/^. The field dependence of the gap 
is in good agreement with the experimental result ob- 
served for copper benzoate.^^ On the basis of quantum 
SG theoryji^ they also provided a good description for 
the resonance field of electron spin resonance (ESR) ob- 
tained by Oshima et al^ 




FIG. 1: Schematic view of the lowest-energy excitations of 
model ([!]) in a nonzero magnetic field for h = (dashed lines) 
and hj^O (solid lines). Gapless excitations at q^O and at 
incommensurate waves q — izqo for h — have finite gaps for 

The elementary excitations in copper benzoate were 
investigated in detail using high-frequency ESR.^^'^^ In 
addition to copper benzoate, the following substances 
are known to be described by the effective model 
given by eq. (P): Yb4As3^--2^ PM-Cu(N03)2-(H20)2 
(PM = pyrimidine)^-^ and CuCl2-2((CD3)2SO),^-^ 
Because the elementary excitations and thermodynamic 
properties of quantum SG systems have been of great 
interest, numerous theoretical investigations have been 
published'^4r45,^ and the theoretical results have been used 
to analyze experimental results on the above-mentioned 
substances. Moreover, new experiments have been pro- 
posed on the basis of theoretical results . ^^i^^" — 

KCuGaFg is an 5* = 1/2 AFHC system, which is rep- 
resented by the effective model given by cq. ([1]) in a mag- 
netic field.— 1^ KCuGaFg is composed of corner-sharing 
CuFg octahedra running along the c axis, as shown in 
Fig.d^i The CuFg octahedra are elongated perpendicu- 
lar to the chain direction, which is parallel to the c axis, 



owing to to the Jahn- Teller effect. Consequently, the hole 
orbitals of Cu^+ ions are linked along the chain direction 
through the p orbitals of F~ ions, which leads to a strong 
antiferromagnetic exchange interaction along the c direc- 
tion of J/fce — 103 K.^^'^a The elongated and compressed 
principal axes of the octahedra alternate along the chain 
direction, as shown by the filled and unfilled bonds in 
Fig.m respectively. This gives rise to the staggered com- 
ponent of the g tensor and the DM interaction with the 
alternating D vector. Because of the c glide plane at 
±6/4, the ac plane component of the Di vector alter- 
nates along the chain direction but the b component does 
not. Thus, the Di vector should be expressed as 

D,^[{-iyD^,Dy,i-iyD,], (4) 

where the x, y and z axes are chosen to be parallel to 
the a* (_L6,c), b and c axes, respectively. For these rea- 
sons, the staggered transverse magnetic field hi is in- 
duced, when this compound is subjected to an external 
magnetic field H. If the y component Dy is negligible, 
then the Di vector is expressed as Di — {~iyD. The 
classical spin arrangements for the staggered and uni- 
form D vectors are the canted Neel state and the spiral 
state, respectively. 




FIG. 2: Crystal structure of KCuGaFe viewed along the a 
axis. Dotted lines denote the chemical unit cell. Filled and 
unfilled bonds respectively denote the elongated and com- 
pressed axes of CuFs octahedra. 

KCuGaFg differs from other quantum SG sub- 
stances in its large exchange interaction of J/fce — 103 K 
and its wide range of the proportionality coefficient, 

= 0.03 -0.174^1^ Thus, KCuGaFe is expected to be 
suitable for obtaining a comprehensive understanding 
of the systems described by model ([1]). In a previous 
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papei)^ we observed as many as about ten ESR modes 
in KCuGaFg and found that most of the ESR modes can 
be very well explained by quantum SG field theory.— 
The magnetic-field-induced gap in KCuGaFg was also 
confirmed by NMR measurement and a gap proportional 
to H^/^ was observed for H\\cr— 

Specific heat measurement is a powerful tool necessary 
for obtaining a comprehensive understanding of the exci- 
tation nature. In the present paper, we report the results 
of specific heat measurements on KCuGaFg in three dif- 
ferent field directions. Experimental results are analyzed 
by performing calculations based on quantum SG field 
theory. As shown below, the specific heat of KCuGaFg 
in a magnetic field can be basically understood in terms 
of quantum SG field theory. However, we found that ad- 
ditional excitation modes, whose origins are not known, 
make significant contributions to the specific heat. This 
paper is organized as follows: in section 2, we summarize 
the theory of elementary excitations and specific heat 
ing the quantum SG model. Experimental details are 
reported in section 3. The experimental results, their 
analyses and a discussion are given in section 4. Section 
5 is devoted to a conclusion. 



II. ELEMENTARY EXCITATIONS AND 
SPECIFIC HEAT IN QUANTUM SG MODEL 
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FIG. 3: Illustration of low-energy excitations of model ((TJ 
around q = G. The soliton, antisoliton, soliton resonance and 
three breathers are labeled as S, S, Es and Mi {i = l, 2 and 
3), respectively. 



The elementary excitations characteristic of quantum 
SG model are solitons, antisolitons and their bound states 
called breathers In field-theoretical language, the ex- 
citation energies are expressed as the masses of these 
quasi-particles. Figure [3] illustrates low-energy excita- 
tions around q = 0. Because of the transverse staggered 
field h induced by the external magnetic field, the gap- 
less excitations at q = 0, tt, ± qo and 7r± qq for h = have 
finite gaps. 

The soliton mass Mg was calculated analytically by 
Essler et aZ.~ as 
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T 



1+e 
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,(5) 



where v is the dimensionless spin velocity, ^ is a param- 
eter given by ^= [2/(7ri?^) — 1]"^ and c is a parameter 
depending on the magnetic field. The field dependences 
of these parameters are given in the literaturer^ilSiS 
v^tt/2, and c^l/2 for H^O. This result is 

applicable in a wide magnetic field range up to the satu- 
ration field given by Hg = 2 J/ gfiB- 

The breathers M„ (n = 1, 2, • • • ) corresponding to the 
excitations at q — and tt have hierarchical structures 
labeled by integer n. The mass of the nth breather is 
given in terms of the soliton mass Mg and parameter ^ 
as 



M„ = 2AfsSin 



(6) 



The number of breathers is limited by n< [^"-"^J.— In a 
low magnetic field of gfinH/ J < 0.5, breathers up to the 
third order can exist. The soliton resonance labeled as 
Eg in Fig. [3] corresponds to the q = excitation on the ex- 
citation branch connected to the soliton and antisoliton. 
Its energy is given by Eg ~ [M^^ 4- {g^iBHfY^'^}'^ 

The specific heat in SG field theory can be obtained 
exactly by taking advantage of its integrability. It has 
been evaluated by two approaches: the thermodynamic 
Bethe ansatz (TBA)— and quantum transfer matrix 
(QTM) methods.-^i^ In the TBA method, a set of in- 
tegral equations is derived from the Bethe ansatz. The 
number of integral equations is finite only at special val- 
ues of the compactification radius, including the SU(2) 
symmetric point R = l/-\/27r. Thus, the application to 
the present problem is practically limited to the SU(2) 
symmetric radius, where the first breather mass Mi de- 
generates with soliton mass Mg. Precisely speaking, how- 
ever, the uniform field H breaks the SU(2) symmetry and 
thus renormalizes the compactification radius R. This 
renormalization effect can be taken into account by us- 
ing the QTM method, which allows the calculation of 
specific heat for an arbitrary radius. These calculations 
of the exact specific heat in SG field theory has been suc- 
cessfully applied the analysis in the experimental data on 
PM-Cu(N03)2-(H20)2 ^^'^^ and Yb4As3,24 which exhibit 
a field-induced gap similar to that in the present system. 
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III. EXPERIMENTAL DETAILS 

KCuGaFg single crystals were grown by the horizontal 
Bridgman method from the melt of an equimolar mixture 
of KF, CuF2 and GaF^ packed into a Pt tube of 10 or 15 
mm inner diameter and 100 mm length. The materials 
were dehydrated by heating in vacuum at about 100°C 
for three days. One end of the Pt tube was welded and 
the other end was tightly folded with pliers. The tem- 
perature at the center of a furnace was set at 750° C, and 
was decreased at a rate of l°C/h to 500° C. Transparent 
light-pink crystals with a maximum size of 10 x 15 x 5 
mm'^ were obtained. These crystals were identified as 
KCuGaFg by X-ray powder diffraction analysis. 

The crystallographic a, b and c axes were determined 
by X-ray single-crystal diffraction. KCuGaFg crystals are 
cleaved along the (1, 1, 0) plane. The magnetic suscepti- 
bilities Xa, Xb and Xc are anisotropic below 50 K because 
of the DM interaction)^ The magnitudes of the suscep- 
tibilities below 50 K decrease in the order Xc> Xb> Xa- 
Thus, the crystallographic axes can be identified from the 
susceptibility measurements. Specific heat measurements 
were carried out down to 0.35 K in a magnetic field of up 
to 9 T using a Physical Property Measurement System 
(Quantum Design, PPMS). Magnetic fields were applied 
along the a, b and c axes. 



IV. RESULTS AND DISCUSSION 

Figure m shows the low-temperature total specific heat 
C'totai measured at zero magnetic field. No magnetic or- 
dering was observed down to 0.35 K, which verifies the 
good one-dimensionality of the present system. Ctotai is 
almost proportional to temperature below 4 K, which is a 
characteristic of the 5 = 1/2 AFHC^''^^^ To evaluate the 
lattice contribution, we used the theoretical specific heat 
of the 5=1/2 AFHC-^i^ The solid line in Fig.S] indi- 
cates the theoretical specific heat Cafhc calculated with 
J/fcB ^ 103 K^^i^ For ksT/J < 0.1, the specific heat of 
the 5 = 1/2 AFHC is approximately expressed as 



Cafhc — 



2RkBT 
3J ' 



(7) 



where R is the gas constant and the spin state is equiv- 
alent to the Tomonaga-Luttinger liquid (TLL). Ctotai is 
the total of the magnetic Cmag and lattice Ciatticc contri- 
butions. The lattice contribution Ciatticc shown by the 
dashed line in Fig. |4] was obtained by subtracting Cafhc 
from Ctotai- The magnetic specific heat in a finite mag- 
netic field was obtained by subtracting Ciatticc from the 
total specific heat Ctotai- 

Figure [S] shows the low-temperature magnetic specific 
heat Cmag obtained at various magnetic fields parallel to 
the c axis. For this field direction, the proportionality 
coefficient Cg = h/H is the largest among the three direc- 
tions and Cs = 0.17. In a finite magnetic field, Cmag in- 
creases exponentially with increasing temperature, which 
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FIG. 4: Low-temperature total specific heat Ctotai measured 
at zero magnetic field. Specific heat data are plotted as open 
circles. Cafhc is the theoretical low-temperature specific heat 
given by eq. (O with J/ks = 103 K, and Ciatticc is the lattice 
contribution estimated by subtracting Cafhc from Ctotai - 



indicates the emergence of a magnetic-field-induced gap. 
With further increasing temperature, Cmag exhibits a 
rounded shoulder and increases linearly. The shoulder 
shifts to a higher temperature and becomes broader as 
the magnetic field increases. This shows that the gap 
increases with the magnetic field. 

We analyzed the experimental results using the ex- 
act specific heat in SG field theory. To verify the accu- 
racy of the calculation and to determine the effect of the 
field renormalization of R, we compared the experimental 
data with the theoretical specific heats obtained by the 
TBA and QTM methods. In the TBA calculation, R is 
fixed to the SU(2) symmetric radius as discussed above. 
In the QTM calculation, R is obtained as a function of H 
using the exact Bethe ansatz solution of the Heisenberg 
chain without the staggered field Thus, the only ad- 
justable fitting parameter is the soliton mass Mg in both 
calculations. 

The solid lines in Fig. [5] show the fits by Ctba(A) with 
the soliton mass shown in Fig. [Bl On the whole, the ex- 
perimental specific heat for H \\c is well described by the 
theory based on quantum SG field theory, although some 
discrepancy is observed. The following three factors may 
have led to the discrepancy between the experimental 
and theoretical results: the breaking of the SU(2) sym- 
metry, the error in the estimation of the lattice contribu- 
tion Ciatticc and the effect of unknown excitation modes 
observed in previous ESR measurements.— 

The effect of breaking the SU(2) symmetry is deter- 
mined by comparing the TBA calculation with R fixed 
at 1/^/2tt and the QTM calculation with the renormal- 
ized R. Figure [7] shows a comparison between the re- 
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FIG. 5: Temperature dependence of magnetic specific heat 
Cmag measured at 7^=1,3,6 and 9 T for || c. Each set 
of data is shifted upward by a multiple of 0.5 J/(K mol). 
Open circles indicate experimental data and solid lines are 
the theoretical specific heat Ctba(A) based on quantum SG 
field theory fitted with the soliton mass A shown in Fig.|6] 



suits of analyses based on the TBA and QTM meth- 
ods. As shown in Figs. [5] and [71 the temperature de- 
pendences of the specific heat and the soliton masses 
obtained from both methods are similar. This means 
that, in the present magnetic field range, the effect of 
the renormalization of R is not important and the sys- 
tem can be well approximated by the SU(2) symmetric 
radius R = l/^/2n. In fact, in the present material, the 
exchange interaction J/fce = 103 K is dominant over the 
magnetic field used in the experiments42i^ In addition, 
the specific heat is not very sensitive to a small change 
in R. 

The soliton masses A for iJ 1 1 c evaluated from the anal- 
yses using the TBA and QTM methods are shown in 
Fig. ini as a function of H^/^. It is clear that A is pro- 
portional to H^/^. This relation can be derived by set- 
ting ^ = 1/3 in eq. ([5]) . Note that in the linear spin wave 
theory, the energy gap exhibits the field dependence of 
AccH^/^. Therefore, the field dependence of the soliton 
mass for H \ \c agrees well with that in quantum SG field 
theory. 

The dotted line in Fig. [S] indicates the soliton mass A 
obtained from previous ESR measurements.— The soli- 
ton mass obtained from the specific heat measurements 
is 1.26 times as large as that obtained from the ESR mea- 
surements. The overestimation of the soliton mass is at- 




FIG. 6: Soliton mass A(=Afs) for iif || c as a function of 
H^^^ obtained from the fits of the theoretical specific heat 
calculated using the TBA and QTM methods. The dot- 
ted line shows the soliton mass obtained from the ESR 
measurements.— 
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FIG. 7: Analyses of the specific heat obtained using the TBA 
and QTM methods based on quantum SG field theory. Circles 
show the magnetic specific heat measured &t H = &T for H \\ c. 
Dashed and solid lines indicate the theoretical specific heat 
calculated using the TBA and QTM methods, respectively. 



tributed to the unknown higher-energy excitation modes 
U2 and U3 observed in ESR measurements.'^'^ In another 
quantum SG system PM-Cu(N03)2-(H20)2, the ratio of 
the gap obtained from specific heat measurements^ to 
that obtained from ESR measurements— is similar to 
that observed for KCuGaFg. While the small disagree- 
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ment in the soliton mass is found, the experimental Cmag 
for II c is weU described within the framework of quan- 
tum SG field theory, as shown in Figs. [5] and [71 




T[K] 



FIG. 8: Magnetic specific heat Cmag measured for H \\a. Each 
set of data is shifted upward by a muhiple of 0.25 J/(K mol). 
Dashed lines denote fits by Ctba(A) with a single soliton 
mass A, while solid lines indicate fits using eq. (|8]) with the 
two gaps Ai and A2 shown in Fig. [9l 

Next, we show in Fig.[5]the low-temperature magnetic 
specific heat Cmag obtained at various magnetic fields for 
H II a. For this field direction, the proportionality coef- 
ficient is the smallest and Cs=0.031i^ From the expo- 
nential increase in Cmag, we see that the gap opens in 
a magnetic field. However, because of the small propor- 
tionality coefficient Cg , the magnitude of the gap is much 
smaller than that for || c at the same magnetic field. 
First, we fitted the theoretical specific heat Ctba(A) to 
the experimental data. The soliton mass A obtained is 
plotted in Fig. [HI as a function of H^^^. As shown by the 
dashed lines in Fig. |S1 some discrepancy was observed 
between the experimental and theoretical results. A sim- 
ilar discrepancy was also observed in copper benzoate 
when an external field was applied along the a" direc- 
tion, for which the gap becomes smallest. '^^ Although A 
is proportional to H^^^, as predicted by quantum field 
theory, its magnitude is 1.45 times as large as that ob- 
tained from the ESR measurements. It is considered that 
this discrepancy for KCuGaFg arises from the higher- 
energy unknown mode Ui and multiple excitation modes 
C„, which cannot be explained within the framework of 
quantum SG field theory?S2, 

The energy of the Ui mode observed for _ff || a is larger 




FIG. 9: Soliton mass A (= M^) for ii" || a as a function of H'^^^ 
obtained from fits of the theoretical specific heat Ctba(A) 
(closed squares). Closed and open circles indicate the gaps 
Ai and A2 obtained by fits using eq. (|8]), respextively. Open 
squares indicate the soliton mass estimated by ESR measure- 
ments. 



than the mass of the first breather and is compara- 
ble to those of the second and third breathers. The 
Cn modes are multiple excitations of the Es and the 
breathers, whose energies are given by Es + Mn- Note 
that these modes were also observed in ESR measure- 
ments on PM-Cu(N03)2-(H20)2'^ These excitations can 
contribute to the magnetic specific heat because they are 
excitations from the ground state. Then, assuming that 
the contribution of these higher-energy excitations is ef- 
fectively represented as Ctba(A2) with the secondary 
gap A2, we express the theoretical specific heat as 

Cmag = \ [CtBa(Ai) + CtBa(A2)] . (8) 

The first term expresses the contribution of the excita- 
tions of the quantum SG model, while the second term 
expresses the contribution of other higher-energy exci- 
tations. The primary gap Ai corresponds to the true 
soliton mass. The solid lines in Fig. |S] show fits using 
eq. ([H . The experimental specific heat in various mag- 
netic fields for 7? II a ia well described by C^g. The two 
gaps, Ai and A2, obtained from the fits are plotted as 
functions of H^/^ in Fig. [Altogether with the soliton mass 
estimated from the ESR measurements 1^ Both energy 
gaps are proportional to i?^/^ for > 3 T. The primary 
gap Ai coincides with the soliton mass estimated from 
the previous ESR measurements. This means that the 
lower excitations for -ff || a can be understood within the 
framework of quantum SG field theory. 

Figure [TUl shows the temperature dependences of the 
total specific heat Ctotai for H\\h and -ff || c, for which 
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FIG. 10: Comparison between total specific heat Ctotai for 
H \\b and H \\c measured at various magnetic fields. Each set 
of data is shifted upward by a multiple of 0.5 J/(mol K). 



FIG. 11: Magnetic specific heat Cmag measured for H \\ b. 
Each set of data is shifted upward by a multiple of 0.2 J/(mol 
K). Dashed lines show fits by the theoretical specific heat 
Ctba(A), while solid lines indicate fits using eq. (|8|l. 



the proportionality coefficients are Cg = 0.16 and 0.17, 
respectively.'^" The values of magnetic specific heat for 
these two field directions were expected to be similar, 
because the proportionality coefficients are close to each 
other. However, the values of specific heats are clearly 
different, as shown in Fig. [TOl It is apparent that the 
magnetic- field-induced gap for H \ \bis much smaller than 
that for II c. A notable feature is that the specific heat 
for _ff II 6 is almost linear in temperature for < 2 T. This 
indicates that the ground state is gapless at low mag- 
netic fields. The low-magnetic-field specific heat dose not 
exhibits an anomaly indicative of 3D ordering down to 
0.35 K, although the ground state appears gapless. These 
observations appear to be consistent with the result of 
previous ESR measurements for H \\b^ where we observed 
an intense unknown mode with resonance frequency 
proportional to H - with i?c =i 2.5 T.^ 

Figure [TT] shows the magnetic specific heat measured 
at various magnetic fields for H\ \ h. Fitting of the the- 
oretical specific heat Ctba(A) to the experimental data 
is not successful, as shown by the dashed lines in Fig. [TTl 
and contrary to the case of H \\ a, we obtain a soli- 
ton mass A considerably smaller than that estimated 
from ESR measurements (see Fig. [T2|) . Soliton mass 
A is not proportional to H'^/^ but can be expressed 
as A = A{if2/3_^^(o)2/3} ^j^j^ i/^(0) = 1.22T. In pre- 
vious ESR measurements, we observed three unknown 
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FIG. 12: Soliton mass A(=Ms) for -ff || 6 as a function of 
li^^^ obtained from the fits of the theoretical specific heat 
Ctba(A) (closed squares). Closed and open circles indicate 
the gaps Ai and A2 obtained by fits using eq. respectively. 
Open squares indicate the soliton mass obtained from ESR 
measurements. 
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FIG. 13: Low-temperature total specific heat Ctotai measured 
at various magnetic fields for -ff || 6. Open circles indicate 
Ctotai at zero field. Closed symbols indicate Ctotai in nonzero 
magnetic fields. 
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FIG. 14: Magnetic field dependence of total specific heat 
Ctotai measured at r = 0.4 K for H \\ b. The quantum phase 
transition where a gapless ground state changes into a gapped 
state is observed. 



modes, [/i, U2 and C/s, whose excitation energies are close 
to or lower than the first breather mass added to in the Ui 
modci^ The unexpectedly small soliton mass obtained 
from the present specific heat measurements should be 
due to these low-energy unknown modes. Then, assum- 
ing that the contribution of these unknown modes is ef- 
fectively expressed in terms of the secondary gap A2, 



which is much smaller than the soliton mass Ai, we de- 
scribe the magnetic specific heat by eq. ([5]), as in the 
case of H\\a. The solid lines in Fig. [11] indicate the 
fits using eq. (jS]) with the two gaps shown in Fig. [T^l 
the experimental value of Cmag for H \\h is effectively 
described by eq. dH). Similar to the soliton gap A ob- 
tained from the fit by Ctba(A), both the primary gap 
Ai and the secondary gap A2 obtained for > 3 T are 
expressed as A^ = A{H'^/^ - Hc{if/^} (i = 1 and 2) with 
ffc(l) = 0.51T and i/c(2) = 1.98T. This value of 77c(2) 
is close to the critical field iJc — 2.5T observed in ESR 
measurements for H \\ 

In Fig. [131 we show the total specific heat Ctotai mea- 
sured below 1 K at various magnetic fields. In this tem- 
perature range, the lattice contribution Ciattico is neg- 
ligible. Ctotai appears to be linear in temperature for 

< 3 T, while for > 5 T, Ctotai exhibits exponential 
temperature dependence. This indicates that the gapless 
ground state changes into a gapped state between 3 and 5 
T. To confirm the magnetic-field- induced quantum phase 
transition, we performed a field scan of specific heat at 
T = 0.4K. The result is shown in Fig.[14l The total spe- 
cific heat Ctotai has a peak &i H = 0.6 T and a shoulder 
at H = 3 T. This behavior was not observed for H\\a 
and H\\c. With further increasing magnetic field, Ctotai 
decreases rapidly and becomes almost zero at H = 5.5 T. 
From these observations, we can deduce that the gapless 
ground state for H\\b becomes a gapped state at a criti- 
cal field He of 3 and 5.5 T. The critical field is larger than 
He — 2.5 T observed in ESR measurements. At present, 
the origins of the low- field gapless state ioi H \\b and the 
unknown ESR modes that make a large contribution to 
the specific heat are not clear. The uniform b axis com- 
ponent of the D vector, which gives rise to a helical spin 
structure or soliton lattice in a magnetic fieldjSi may be 
responsible for these effects. 



V. CONCLUSION 

We have presented results for the specific heat of 
KCuGaFg in magnetic fields. It was clearly observed 
that an excitation gap opens in a magnetic field and 
increases with increasing magnetic field. Specific heat 
data were analyzed using two theoretical values of spe- 
cific heat calculated by the TBA and QTM methods, 
both of which are based on the quantum SG model. In 
the TBA method, the SU(2) symmetry is assumed, while 
in the QTM method, the breaking of the SU(2) symme- 
try owing to the external magnetic field is taken into ac- 
count. The specific heat for _ff || c is well reproduced by 
these two calculations with the soliton mass close to that 
obtained from previous ESR measurements. The soli- 
ton mass was found to be proportional to H^^^. These 
results indicate that the thermodynamic properties for 
H \\ c are well described by quantum SG field theory, and 
that the breaking of the SU(2) symmetry is negligible 
in our experimental magnetic field range because of the 
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large exchange interaction of J/fce = 103K. 

In the cases of H\\a and H \\h, we found a signifi- 
cant contribution of additional excitations that cannot 
explained in the framework of quantum SG field theory. 
For these field directions, we analyzed the specific heat 
data using eq. ([8|), in which the contribution of the un- 
known excitations is effectively taken into account by as- 
suming a secondary gap A2. For H \ \a, the primary gap 
Ai agrees with the soliton mass Aesr obtained from the 
ESR measurements, while for H \\ b, Ai does not agree 
with Aesr- We found that for H\\ b, the ground state 
is gapless at low magnetic fields, and that with increas- 
ing magnetic field, a quantum phase transition occurs 
between gapless and gapped ground states. 
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